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Abstract:Neutrinos and antineutrinos of all three flavours are emitted during the post 
bounce phase of a core collapse supernova with Vu,/v t (Pu,/vt) having average energies more 
than that of v e {p e ). They can be detected by the new earth bound detector like SNO and 
Super-Kamiokande which are sensitive to neutrinos of all three flavours. In this letter we 
consider the effect of flavour oscillations on the neutrino flux and their expected number 
of events at the detector. We do a three-generation analysis and for the mass and mixing 
schemes we first consider the threefold maximal mixing model consistent with the solar and 
the atmospheric neutrino data and next a scenario with one Am? ~ 10 _11 eV 2 (solar range) 
and the other Am? ~ 10~ 18 eV 2 , for which the oscillation length is of the order of the super- 
nova distance. In both these scenarios there are no matter effects in the resultant neutrino 
spectrum and one is concerned with vacuum oscillations. We find that though neutrino os- 
cillations result in a depletion in the number of v e and v e coming from the supernova, the 
actual signals at the detectors are appreciable enhanced. 
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The question whether neutrinos are massive or not has been answered. After 537 
days of data on atmospheric neutrinos the Super-Kamiokande has finally confirmed the 
existence of non-zero oscillations and hence mass for the muon neutrinos Q . The Super- 
Kamiokande (SK) data confirmed the depletion in the atmospheric muon neutrino flux 
which the Kamiokande, 1MB and the Soudan experiments had observed before. At 
90% CL. the mixing parameters allowed by all the atmospheric neutrino experiments 
combined are @ 5x 10~ 4 eV 2 < Am 2 < 6x lO~~ 3 eV 2 , sin 2 26 > 0.8 for the v^—v T channel 
and 10" 3 < Am 2 < 7 x 10~ 3 eV 2 , sin 2 26 > 0.8 for the - v s (Am 2 > 0) oscillation 
mode, while the oscillation parameter region for the — v e mode is completely ruled out 
by the data from the CHOOZ experiment [Q. The other puzzle which has warranted 
neutrino oscillations as a possible solution is the solar neutrino deficit problem. The 
three solar neutrino detectors, the Homestake, Kamiokande and Gallex (also Sage) 
have been observing neutrino flux far less than that predicted by the standard solar 
model |4]]. This deficit can be explained by neutrino oscillations in vacuum for 5m 2 ~ 
0.615 x lO -10 ^ 2 and sin 2 26 ~ 0.864 [BJ or by MSW resonant flavour conversions 
U for Am 2 ~ 5.4 x 10~ e eV 2 and sin 2 26 ~ 7.9 x 10~ 3 (non-adiabatic solution) and 
Am 2 ~ 1.7 x 10~ 5 eV 2 and sin 2 29 ~ 0.69(large angle solution) 0. The first results from 
the SK solar neutrino flux measurements favor the long wavelength vacuum oscillation 
solution with large angle mixing ||. When the SK solar v data is combined with the 
earlier data then the corresponding values are Am 2 ~ 6.5 x 10~ ll eV 2 , sin 2 29 ~ 0.75 
for vacuum oscillations and Am 2 ~ 5 x 10~®eV 2 , sin 2 29 ~ 5.5 x 10~ 3 for non-adiabatic 
MSW resonant flavour conversion 0. The threefold maximal mixing model JlO], |ll 



can explain both the solar and the atmospheric neutrino data simultaneously |12j. In 
a recent paper [jnj it has been shown that the maximal mixing model can account for 
both the SK atmospheric data and the CHOOZ data provided the relevant Am 2 is in 
the range 4 x l0~ 4 eV 2 < Am 2 < 1.5 x 10~ 3 eV 2 , along with the solar v data. One hopes 
to find a definite solution to the solar neutrino problem once the Sudbury Neutrino 
observatory (SNO) which is the first heavy water detector becomes operational ]TJ . 
Though one of the principal motivation for these two detectors was to throw light on 
the solar neutrino problem but they are equally useful for detecting the neutrinos from 
a nearby supernova event. 
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The core of a massive star (M > SM & ) starts collapsing once it runs out of nuclear 
fuel. The collapse continues to densities beyond the nuclear matter density after which 
a bouncing of the in falling matter takes place leading to supernova explosion and 
the formation of a protoneutron star. Only a small fraction of the huge gravitational 
energy released in the process goes into the explosion and all the rest of the energy 
is carried away by neutrinos and antineutrinos of all three flavours. These neutrinos 
for galactic supernova events can be detected by detectors like the SNO and SK. In 
contrast to the solar, the atmospheric as well as the accelerator/ reactor neutrinos 
where one has neutrino flux of a single flavour at the source, postbounce supernova 
neutrinos (antineutrinos) start from the source in all three flavours but with v^jv T 
[p^j u T ) having average energies more than that of z/ e (z/ e ) and it is an interesting problem 
to study whether their flux and their signal at the terrestrial v detectors get appreciably 
altered in reaching the earth if neutrinos do oscillate. In this work we give quantitative 
predictions for the number of neutrino events coming from a typical type II supernova 
at a distance of lOkpc in both SNO and SK and show how the number of events for 
each detection process would change in case oscillations do take place. 

There have been various attempts before to estimate the effect of non-zero neutrino 
mass and mixing on the expected neutrino signal from a galactic supernova. Matter 
enhanced resonant flavour conversion has been observed to have a large effect on the 



v e signal |J| [16|, ^7|. The v e events of course remain unchanged in this case. With 



vacuum oscillations we can expect an increase in both the v e and v e signal. Burrows 



et al. |16| have considered for SNO, the effect of vacuum oscillations as well and have 
found that with two-flavours the effect of vacuum oscillations on the signal is small, 
using their model predictions for the different v luminosities. 

We have considered a three-generation mixing scheme and have calculated the 
effect of neutrino oscillations on the signal from a 20 M Q supernova model devel- 



oped recently [19]. First we do our calculations for the threefold maximal mixing 
model consistent with the solar (Am 2 ~ 10~ u eV 2 ) and the atmospheric neutrino data 
(Am 2 ~ 10~ 3 eV 2 ). For Am 2 ~ !0~~ 3 eV 2 normally we expect matter enhanced reso- 
nance in the supernova. But for the particular case of maximal mixing it has been 
shown before, both numerically |2(J and analytically |21 |, that there are absolutely no 
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matter effects in the resultant neutrino spectrum on earth. Though the arguments in 
both these previous papers are for solar neutrinos, extension to the case of supernova 
neutrinos is straightforward. Hence in this scheme we are concerned with vacuum os- 
cillations. We also consider a second scenario where we take one of the mass square 
differences in the solar vacuum oscillation solution range while the other is ~ lO~ 18 eV 2 
where the oscillations wavelength A ~ L, the distance of the supernova from the earth 
and hence oscillations are observable in the neutrino spectrum. In this scenario of 
course there is no chance of a MSW resonance and we have vacuum oscillations. We 
find appreciable enhancement in the expected u e and v e charge current events for both 
SNO and SK in both scenarios even with vacuum oscillations. 

The differential number of neutrino events at the detector for a given reaction 
process is 



One uses for the number of neutrinos produced at the source N v (t) = L v (t) / {E v (t)) 
where L u (t) is the neutrino luminosity and (E u (t)) is the average energy. In (H) a(E) 
is the reaction cross-section for the neutrino with the target particle, L is the distance 
of the neutrino source from the detector(lOkpc), n is the number of detector particles 
for the reaction considered and f u {E) is the energy spectrum for the neutrino species 
involved. For the neutrino luminosity and average energy we use the values of Totani 
et al. |19| for a 20 M & type II supernova model based on the hydrodynamic code 
developed by Wilson and Mayle. Though in their paper Totani et al. observe that 
the neutrino spectrum is not a pure black body, but we as a first approximation use 
a Fermi-Dirac spectrum for the neutrinos, charaterised by the v temperature alone 
for simplicity. The effect of a chemical potential is to cut the high energy tail of 
the neutrino spectrum and we also study it's effect on the the v signal and on the 
enhancement of the signal when oscillations are introduced. We find the v signal for 
the various detection processes as a function of energy by integrating out time from 
([I]). By integrating over energy as well we get the total number of events for the 
reaction concerned. These are the expected number of events. 

In the presence of oscillations of massive neutrinos more energetic v^ipy) and ^ T (^r) 
get transformed into u e (u e ) which modifies the numbers that we obtain using (||). The 
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general expression for the probability that an initial v a gets converted to a Up after 
traveling a distance L in vacuum is 

tyL 

Pu a u = 5 a (3-AY J U ai U pi U aj U Pj sin 2 — (2) 

j>i ^ij 

where a = e, /z, r, .. and i, j = 1, 2, 3, .. 
• A i7 = 2.5 x 10~ 3 A;m J> v fl 



Am^j = m 2 — mf 



U a i are the components of the mixing matrix. For the mass and mixing parameters we 
consider two scenarios. 

• scenario 1: Here we consider threefold maximally mixed neutrinos with the 
mass spectrum Am^ ~ Am^g ~ lQ~ z eV 2 corresponding to the atmospheric 
range while Am^ ~ 10~ ll eV 2 in accordance with the solar neutrino problem. 
The oscillations due to all the mass differences are averaged out to 1/2 as A << L, 
and hence the expression for the various probabilities in this case relevant for us 
are 



11. 12 



Pu e u e — g (3) 

PuuU e 4" Pv T u e 1 Pv e v e (4) 



We call this Case 1. 



scenario 2: Here we set Amf 2 ~ 10~ 18 eV 2 for which A ~ L and the oscillation 
effects are observable while Am 2 3 « Am^ ~ 10~ u eV 2 (solar range). If we con- 
sider the Maiani parametrisation of the mixing matrix U |22| then the expression 
for the probabilities are 

tyL 1 

Pv„v e = 1 - sin 2 26i 2 cos i 9i 3 sin 2 - -sin 2 26 13 (5) 

A12 2 

-PuuVe 4" Py T U e 1 Pu e U e (6) 
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For this case the oscillations due to Am 2 3 and Am| 3 are averaged out as the 
neutrinos travel to earth but those due to Am 2 2 survive. For 9 is we consider 
two sets of values allowed by the solar v data. We have done our calculations 
for sin 2 2# 13 = 1.0 (the maximum allowed value) and with sin 2 2# 13 = 0.75 (the 
best fit value) ||. The first set is called Case 2a while the second is called Case 
2b. Since nothing constrains Am 2 2 in this scenario we can vary 9ii and study 
it's effect on the v signal. We have tabulated our results for sin 2 2# 12 = 1.0 since 
it gives the maximum increase in the signal from the no oscillation value. 

The corresponding expressions for the antineutrinos will be identical. We note that 
because the energy spectra of the z/ M and v T are identical, we do not need to distinguish 
them and keep the combination P v Ve + P VtV& . We have made here a three-generation 
analysis where all the three neutrino flavours are active. Hence if both the solar v 
problem and the atmospheric v anomaly require v oscillation solutions, then in the 
scenario 2, the atmospheric data has to be reproduced by — v s oscillations. We 
are interested in this scenario as only with neutrinos from a supernova can one probe 
very small mass square differences ~ 10~ 18 eV 2 . To find the number of events with 
oscillations we will have to fold the expression ([!]) with the expressions for survival and 
transition probabilities for the neutrinos for all the cases considered. 

In Table 1 we report the calculated number of expected events for the main reactions 
in H 2 and D 2 0. Column 2 of Table 1 gives the expected numbers for the model under 
consideration when the neutrino masses are assumed to be zero. Column 3,4,5 give the 
corresponding numbers for the two neutrino mixing scenarios that we have considered 
(see Table 1 for details). All the numbers tabulated have been calculated for 1 kton of 
detector mass. To get the actual numbers we have to multiply these numbers with the 
relevant fiducial mass of the detector. The efficiency of both the detectors (SNO and 



SK) is taken to be 1 [19|, [23], |24|. The energy threshold is taken to be 5 MeV for both 
SK |23j and SNO |2Hfl . For the cross-section of the (v e — d\ [p e — d), {v x — d) and (v e —p) 
reactions we refer to [53]. The cross-section of the {y e (pe) — e ~) and (v x — e~) scattering 
has been taken from |27j while the neutral current (v x — 16 O) scattering cross- section 



is taken from [ 23[ . For the 16 0(z/ e — e ) W F and 16 (^ e , e + ) 16 iV reactions we refer to [26 



where we have used the cross-sections for the detector with perfect efficiency From a 
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comparison of the predicted numbers in Table 1, it is evident that neutrino oscillations 
play a significant role in supernova neutrino detection. For the neutral current sector 
the number of events remain unchanged as the interaction is flavour blind. 

The 32 kton of pure water in SK detects neutrinos primarily through the capture 
of v e on protons (u e p — > ne + ) and (v e {p e ) — e~) scattering. The energy threshold for 
16 0(z/ e , e~) 16 F is 15.4 MeV and that for 16 0(z/ e , e + ) 16 N is 11.4 MeV, hence these reac- 
tions are important only for very high energy neutrinos. The typical average energies 
of v e and v e from a type II supernova is about 11 MeV and 16 MeV respectively, so we 
do not expect significant contribution from these two reactions. This is evident from 
Table 1 where the 16 events are only 2.1% of the total charge current signal at SK. As 
a result of mixing the mu and tau neutrinos and antineutrinos oscillate (with average 
energy ~ 25 MeV) into v e and v e during their flight from the galactic supernova to 
the detector resulting in higher energy v e and v e and the number of w O events are 
increased appreciably (for Case 1 iy e — 16 O) events go up by 13 times) so that after 
oscillations they are 7% (Case 1) of the total charge current events at SK. The effect 
of oscillations on the (^ e -p) capture is to enhance the expected signal by about 25% 
(Case 1). In all previous studies where the effect of MSW transition on the neutrino 



signal has been studied JT5, [IT], there is no enhancement in the number of expected 
events for the (^ e -p) sector while we do get a significant change in the expected signal 
with vacuum oscillations. For the (z-'e(^e) — e~) scattering the effect of oscillation is 
very small. 

The SNO is the world's first heavy water detector made of 1 kton of pure D2O 
surrounded by ultra pure H 2 0. There are 10 4 phototubes around this entire volume 
which can view only the inner 1.4 kton of water efficiently P3|. We find about 99% 
increase in [y e — d) events and about 46% increase in {p e — d) events for the Case 1. 
From the column 2 of Table 1 we can see that there are more [p e — d) than (v e — d) 
events even though there are more v e than v e coming from the supernova. This is 
because the reaction cross-section a ~ E 2 ' 3 and the v e spectrum is harder than the v e 
spectrum. This also results in a greater enhancement due to oscillations for the {y e — d) 
events, as the difference between the energies of the v e and z/^(z/ T ) is greater than those 
between v e and v^iv^) and hence the effect on the v e events is more. As a result after 
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oscillations are switched on the number of (v e — d) events supersede the (P e — d) events. 
We observe a similar effect for the 16 events, where the v e signal without oscillations 
is more than the v e signal, while the effect of oscillations is more for the latter. The 
effect is more magnified in this case due to the very strong energy dependence of the 
reaction cross-section and also due to the fact that the energy threshold for (P e — 16 O) 
event is lower than for the [y e — 16 O) event. In Fig. 1 we plot the signal due to the 
{v e — d) events as a function of energy , without oscillations and with oscillations for 
the Case I and Case 2b. All the features mentioned are clearly seen. The plot for the 
Case 2b clearly shows oscillations. 

In Fig. 2 we plot the cumulative fluence of the v e coming from the supernova at 10 
kpc without oscillations and with oscillations for Case 1 and Case 2b. It is seen that the 
result of oscillation in fact is to reduce the total number of v e . Yet as seen from Table 
1, we have obtained significant increase in the (v e — d) events and the (v e — 16 O) events. 
The solution to this apparent anomaly lies in the fact that the cross-section of these 
reactions are strongly energy dependent. As a result of oscillations the v e flux though 
depleted in number, gets enriched in high energy neutrinos. It is these higher energy 
neutrinos which enhance the v signal at the detector. This also explains the difference 
in the degree of enhancement for the different processes. For the (u e — d) and (y e — 16 O) 
events, especially for the latter, the effect is huge while for the {y e — e~) scattering it is 
negligible as it's reaction cross-section is only linearly proportional to E. Due to their 
high energy dependent a the 16 0(z/ e , e~) w F events turn out to be extremely sensitive 
to oscillations. A similar argument holds true for the case of the antineutrinos, only 
here the effect of oscillations is less than in the case for the neutrinos as the difference 
between the energies of the v e and v ii jv T is comparatively less as discussed earlier. 

For the scenario 2 we have studied the effect of the mixing angles on the signal. 
For a fixed #13 the effect of oscillations is enhanced if we raise Q^i- The effect of #13 
is more subtle. The effect of oscillations increase with # 13 initially and then decrease. 
We have also checked the effect of a chemical potential fi on the neutrino signal. A 
non-zero \x cuts the high energy tail of the neutrino signal as a result of which the total 
signal goes down for both with and without oscillations, the effect being greater for 
the more energy sensitive reactions. 
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With the supernova model of Totani et al. ]T9|], we have obtained oscillation effects 



in the expected v signal which are significantly larger than those obtained by Burrows 
et al. In the model that Burrows et al. use in their study, the v luminosities L v 

are more than those for Totani et al. model, but the average energy is much smaller, 
particularly for the v e and ^, t (^,t)- Hence their spectra lacks in high energy 
neutrinos which results in almost negligible effect of oscillations in their case. Again 
in the model of Burrows et al. the average energies decrease with time while in the 
model of Totani et al. not only the average energies but also the difference between 
the average energies of v e [pe) and z/^ iT (z/^ iT ) increases with time. The effect of all these 
is to magnify the effect of oscillations in our case. 

In conclusion, we have shown that with the model of Totani et al. even with vacuum 
oscillations we obtain appreciable enhancement in the expected v signal in SNO and 
SK even though the number of neutrinos arriving at the detector from the supernova 
goes down. In contrast to the case where we have MSW resonance in the supernova, 
with vacuum oscillations we get enhancement for both v e as well as v e events. If we 
have a galactic supernova event in the near future and if we get a distortion in the 
neutrino spectrum and an enhancement in the signal, for both v e as well as v e then 
that would indicate vacuum neutrino oscillations. 
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Table 1 The expected number of neutrino events for a 1 kton water cerenkov detector 



reaction 


signal 
without 
oscillation 


signal with oscillation 


scenario 1 


scenario 2 


Case 1 


Case2a 


Case2b 


u e + d^p + p + e~ 


78 


155 


150 


153 


v e + d-^n + n + e + 


93 


136 


133 


135 


v x + d — > n + p + v x 


455 


455 


455 


455 


v e + p — > n + e + 


263 


330 


326 


329 


v e + e~ — > u e + e _ 


4.68 


5.68 


5.61 


5.66 


P e + e~ — > z/ e + e _ 


1.54 


1.77 


1.76 


1.77 


v^rip^r) + e~ -> v^ T {v^ T ) + e" 


3.87 


3.55 


3.50 


3.53 


z/ e + 16 -> e~ + 16 F 


1.13 


14.58 


13.78 


14.45 


P e + 16 -> e+ + 16 AT 


4.57 


10.62 


10.23 


10.53 


+ 16 -> ^ + 7 + * 


13.6 


13.6 


13.6 


13.6 
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Figure Captions 



Fig. 1 The [y e — d) signal at SNO vs neutrino energy without and with oscillations 
for the Case 1 and Case 2b 

Fig. 2 The cumulative v e fluence as a function of the neutrino energy without 
and with oscillations for the Case 1 and Case 2b. Also shown is the fluence for 
comparison. 
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